Azimuthal asymmetries in unpolarized Drell-Yan events on nuclear targets 
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We show that for Drell-Yan events by unpolarized hadronic projectiles and nuclear targets, az- 
imuthal asymmetries can arise from the nuclear distortion of the hadronic projectile wave function, 
typically a spin-orbit effect occurring on the nuclear surface. The asymmetry depends on quantities 
that enter also the spin asymmetry in the corresponding Drell-Yan event on polarized free nucleonic 
targets. Hence, this study can be of help in exploring the spin structure of the nucleon, in particular 
the transverse spin distribution of partons inside the proton. All arguments can be extended also to 
antinucleon projectiles and, consequently, apply to possible future measurements involving nuclear 
targets at the foreseen HESR ring at GSI. 
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I. INTRODUCTION 

High-energy collisions of hadrons can be a very interesting source of information for testing QCD in the nonpertur- 
bative regime. In fact, in the past years several experiments have been performed that still await for a correct inter- 
pretation of the resulting data. Among others, hadronic collisions on transversely polarized proton targets 
where an azimuthally asymmetric distribution of final-state products (with respect to the normal of the production 
plane) is observed when flipping the transverse spin of the target or of the final products, the so-called transverse spin 
asymmetry. Perturbative QCD cannot accommodate for such asymmetries, sometimes as large as 40% also at high 
energy 0. This problem has triggered a growing interest in this field of hadronic spin physics; new data and rapid 
developments are foreseen in the near future, as it emerged in a recently devoted workshop (for a review covering 
also processes with lepton beams, see Refs. 0,0). 

Here, we will concentrate on a series of measurements involving nuclear targets, namely high-energy collisions of 
pions and antiprotons on various unpolarized nuclei 0, |^ ^| where the cross section for Drell-Yan events shows 
an unexpected largely asymmetric azimuthal distribution of the final lepton pair with respect to the production plane. 
It is usually convenient to study the problem in the so-called CoUins-Soper frame of Fig. A lepton plane is 
defined by the direction of emission of the lepton pair and by the z axis which is approximately formed by the average 
of the directions of the colliding hadron momenta Pi and P2 (for a more rigorous definition, see Ref. ^3)- Azimuthal 
angles are measured in a plane perpendicular to z and containing the axis h = q.j,/|q^|, where is the transverse 
momentum of the final lepton pair detected in the solid angle {9, (p). 




FIG. 1: The CoUins-Soper frame. 
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At leading order in as , the cross section |^ 0, ^| can be parametrized as 

= — f 1 + A cos'^ + a sin^ 9 cos (h -\ — sin^ 9 cos 2d> 




(1) 



with A ~ 1 and /i ^ ^ 30 %. Both Leading Order (LO) and Next-to-Leading Order (NLO) perturbative QCD 
calculations give A 1 and /i ^ ~ , because they are based on the assumption of massless quarks and coUinear 
annihilation. In fact, this result holds only if the annihilation direction is well defined and coincides with the z axis of 
Fig.n because the cross section would depend only on B. In general, the assumption of massless quarks leads to the so- 
called Lam- Tung sum rule \ — \ = 2v ^] , which is badly violated by data. More complicated mechanisms, like hi ghe r 
twists or factorization breaking terms at NLO, are not able to explain the size of in a consistent picture pii llTl [l8l| . 
A promising interpretation relies on the observation that at leading twist the cross section for unpolarized Drell-Yan, 
when kept differential also in dq^, contains a term proportional to cos2^(ft,^ ® ^ specific convolution of 

the distribution function that describes the transverse momentum distribution of transversely polarized partons 
inside unpolarized hadrons 'l3l|, while describes the annihilating antiparton partner. 

However, it is important to stress that the above results are obtained using nuclear target: before getting to the 
collision with a bound nucleon that generates the Drell-Yan event, the hadronic beam crosses the nuclear surface 
and travels inside the nuclear medium for some length depending on its energy. The aim of this paper is to try to 
understand whether effects due to the beam-nucleus interaction may assume a relevant role in contributing to the 
observed azimuthal cos 20 asymmetry. In other words, if an intrinsic transverse momentum distribution can affect the 
elementary annihilation of partons leading to an asymmetric azimuthal distribution of final lepton pairs through the 
function , then a hadron deflected by nontrivial Initial-State Interactions (ISI) with the target nucleus can transfer 
its transverse momentum to the final lepton pair and contribute to the same azimuthal asymmetry. More, we will see 
that such ISI can originate azimuthal asymmetries also when no asymmetry is produced at the level of a free proton 
target. 

The highly energetic pions and antiprotons used in measurements reported in Refs. l5lfll lT(llTll | demand a description 
of their internal structure in terms of elementary degrees of freedom (quarks and giuons) in order to correctly describe 
their propagation and to study effects that can modify high-energy azimuthal asymmetries at RHIC or LHC (see, 
for example, Refs. 0,|23|)- Here, we will consider Drell-Yan processes at lower energies of interest for the foreseen 
HESR ring at GSI, where collisions of (polarized) protons and antiprotons with energy in the range 10 30 GeV will 
be studied [2llf2^ I2M 123. 125l l2^ l2^ . Hence, we will adopt a description based on "traditional" degrees of freedom 
as nucleons and pions; this assumption is justified in the following Sec. ^ 

The paper is organized as follows. In Sec. the nuclear damping of the projectile wave function, based on 
the eikonal approximation, is reviewed and the induced modifications to the scattering amplitude are discussed. In 
Sec. mil the spin-orbit effect on the nuclear surface is discussed and the formulae are implemented with an additional 
damping that depends on the polarization component of the incoming hadron beam. In Sec. IIVI it is shown how all 
this can produce azimuthal asymmetries in Drell-Yan events from unpolarized hadron-nucleus collisions even when 
no asymmetry is generated using a free proton target. Finally, in Sec. the discussion is presented about the issues 
of the size of this nuclear effect and of the possibility of disentangling it from genuine asymmetries related to the 
Drell-Yan hard event. Results are summarized in Sec. IVII 



As already stressed in Sec. U we will consider Drell-Yan processes involving nuclear targets at energies of interest 
for the foreseen HESR ring at GSI. At present, the discussion is still open about the two options of an antiproton 
beam of energy 15 GeV hitting a fixed proton target or colliding on a proton beam of energy around 3 GeV, such as 
the center-of-mass energy ranges from 5 to approximately 14 GeV |2lL 122. l2a. xA l25l l2a. l27|. At this scale, the use 
of traditional "nuclear" degrees of freedom (nucleons and pions) is still appropriate for describing the propagation of 
a hadron inside the nucleus, also because spin-orbit effects are still significant. Moreover, a distance in space (the 
so-called formation length) must separate the Drell-Yan hard event from the ISI rescatterings in such a way that the 
two processes do not influence each other; the formation length increases with energy but it is already significant at 
15 GeV beams. In fact, the same problem occurs in the semi-inclusive quasi-elastic electroproduction of protons on 
nuclei about the Final-State Interactions (FSI) of the proton travelling inside the residual nuclear medium after the 
hard event of virtual photon absorption. The NE18 experiment about the ""^^0(6, e''p)X reaction has shown that, 
at 10 GeV energies, the ejectiles interact with the nuclear medium exactly as "normal" protons, before emerging and 
being actually detected as protons. 

At the same time, a projectile with energy 15 GeV impinging on a target nucleus is energetic enough such that the 
collision can be safely considered in the Born Approximation. The main effect of the nuclear medium is, therefore. 



II. SPIN-INDEPENDENT NUCLEAR DAMPING 



to distort the projectile wave function. The most popular and widely adopted approach to such a situation is the 
Glauber method [22, which has a long well established tradition of successfuU results in the field of high-energy 
proton-nucleus elastic scattering (for a review, see also Ref. i^flij)- Again, a parallel with quasi-elastic A{e,e'p)X 
reactions can be established , si nce above the inelastic production threshold, i.e. for proton momenta p 1 GeV/c, 
it has been shown [sil Is^ 13^ that this method gives the same FSI effects as the usual approach based on the 
Distorted Wave Impulse Approximation (DWIA), where the projectile wave function is expanded in partial waves and 
a Schrodinger equation is solved for each partial wave L with a complex optical potential up to a Lmax depending 
on the projectile energy and satisfying a stringent convergency criterion (for a review, see Ref. js^l)- Actually, it 
was directly checked that the equivalence of these methods holds when the Glauber scattering wave is replaced by 
an eikonal wave function ISfij l . Because of the similarity between the FSI problem in A{e,e'p) reactions and the ISI 
problem in p — ^ collisions, we will adopt here the same eikonal approximation for the projectile wave function. 

Since for a fastly moving object the nuclear density can be considered roughly constant but for a small portion 
corresponding to the nuclear surface, the eikonal wave function of the projectile with momentum p can be represented 
by the damped plane wave 

iP{r) K e-P'-'^e'P '' e-P'-" , (2) 

which corresponds to an approximate solution of the Schrodinger equation inside homogeneous nuclear matter with 
|p/| ~ Pi proportional to the ratio between the absorptive (imaginary) potential and |p| = p; it is equivalent to neglect 
the solution of the Schrodinger equation that propagates backward in space (29j. For a non- homogeneous nucleus, the 
same approximation leads to the Glauber wave function. The factor exp(— p/ • R) is due to the incoming proton wave 
being properly normalized to 1, where R is a constant vector whose modulus gives the nuclear radius. A reasonable 
value for pi can be estimated by reproducing, among other things, the NE18 data for the ^^C(e, e'p) reaction 
where the flux of the outgoing proton with p — lA GeV/c is overall reduced by 40%, including the behaviour 
of the related transparency coefficient [s^. It turns out that pi ^ p, for the above NE18 kinematics typically 
Pj ~ 50 MeV/c. Nevertheless, this damping is able to generate "geometric" FSI responsible for the asymmetric 
proton knockout with respect to the scattering plane in exclusive ^'^C{e,e'p) reactions with a polarized lepton beam 
(the so-called fifth structure function problem) |35| . 

The eikonal approximation of Eq. (jSJ is, by construction, a high-energy approximation and its reliability increases 
with the projectile energy, ideally corresponding to the limit of an infinite number of partial waves for the solution in 
the DWIA. In these conditions, as a first step we neglect spin constributions to Eq. (|2Jl (but we will reconsider the 
problem later). Moreover, if the projectile propagates along the z direction, we will assume that also the damping is 
uniformly concentrated along the same direction, i.e. p/ || p = pz; consequently, 

V'(r) « e-P'" e'P' e-P"' = tlj{z) . (3) 

This introduces a (small) error in the asymptotic angular distribution of scattered projectiles, but to our purpose of 
considering only hard events inside the target nucleus it is irrelevant. 

The function (j2Jl can be considered as a plane wave with the complex momentum p-f ip/ = (p+ipi) z — P z = P |35l |. 
Therefore, as a simple plane wave gets Fourier transformed to a (5 distribution, similarly the "complex" plane wave 
exp (iPz) can be Fourier transformed by extending the definition of the 6 distribution to the complex plane. Following 
Ref. [s^Ij since the S distribution of a real argument can be defined as 

S{x — x) = — lim — = lim -— ( z : : ) , (4) 

TT e^O [X — xY -f e^O Im \X — X — It X — X + It J 

the analytic continuation to a complex argument becomes 

5{z -z)= 6{z -{x + ly)) = ^ ( ^ ^ ^ = ^ i — ■ (5) 

2m \z — X — ly z — x + ly J 2m \z — z z — z* J 

In other words, when a function f{x) of the real argument x is convoluted with 5{x — x), the integrand is pinched 
between the two singularities x — x:Lit, which eventually collapse to the same value in the limit e ^ 0. For a function 
/(z) of a complex argument z, the effect is similar if f{z) is analytic and f{z) —>■ for |z| — > oo. In fact, in this case 
we have 



dz6{z-z)f{z)^ f{z) , (6) 



c 



where C is an integration contour in the complex plane extending to Re(z) ±oo on the real axis and closing in the 
upper plane for Im(z) = y > (and viceversa). 



How to interpret the S distribution with a complex argument? With no damping, the eikonal wave in momentum 
space is represented by S{k — p), a singular distribution in k peaked on the momentum p of the projectile. With 
a nuclear damping [pi ^ 0), the momentum distribution is still peaked around p but is no longer singular, since 
the potential inside the nuclear volume distorts the wave function and introduces a width Ap ^ 1 /R around p. For 
energies above the inelastic proton-nucleon threshold, the nuclear potential is mainly imaginary and generates an 
absorption of the projectile flux (here represented by the damping p/); hence, Ap ^ where z is the absorption 
length if the projectile is propagating along the z axis. If z < i?, then Ap is mostly given by the damping, i.e. by pj. 
In conclusion, the distribution 5{k — p — pi) is peaked around p but has a finite width Ap, or equivalently it describes 
a projectile with a finite absorption length z ~ 1/Ap, that is mainly given by the damping p/. 

In order to understand the modifications inside the scattering amplitude, it is useful to consider first the case 
without damping; for the projectile-nucleus collision we have 

a(p) = j dqdk^/';(q) A(q,k) « j dqdk(5(q -pi) yl(q,k) = j dkA(pz,k) , (7) 

where ■(/'p represents the incoming wave function in momentum space centered around p, the function A contains all 
the dynamics of the hard collision, and the integration upon dk covers the portion of phase space complementary to 
the one pertaining the incoming projectile, i.e. it represents the integration upon momenta of nucleons participating 
to the collision and of final lepton products. Then, if the function A is well behaved and vanishes for aymptotically 
large momenta, Eq. lO can be generalized to 

a(p-t-ip/) = J dQdkS{q~P)A{Q,k) ^ J dkA{{p + ipi)z,k) . (8) 
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FIG. 2: The leading-twist contribution to the Drell-Yan process. 



If the invariant mass M of the Drell-Yan lepton pair is kinematically constrained in a range where the annihilation 
proceeds through a virtual photon, the dominant contribution at leading twist is represented by the handbag diagram 
of Fig. 121 [13 , where the blobs represent all the possible residual states of the two hadrons after the hard event that 
are not observed. Hence, a sum over all possible configurations is understood in the two cuts. As such, the process 
is by construction semi-inclusive and the cross section will be proportional to an incoherent sum of several transition 
probabilities, each one leading to the same final lepton pair. Neglecting the irrelevant phase space coefficient, the 
cross section with no nuclear damping will be related to 



(9) 



where the number of terms depends on the model used to describe the blobs in Fig. El 

From Eq. JHJ, it is evident that the damping of the projectile wave function implies a shift in the argument of the 
elementary scattering amplitude, namely a„(p) —>■ an{p + ipi). Since p/ p, we have 



a„(P = p + ipi) = a„(p) + ipi^^ 

oF 



Pi=Q 



+ o(p/) , 



(10) 



so that Eq. © gets modified into 



W{p,pi) = W°{p) - 2p/Im[G(p)] -f o{pi) , 



(11) 



with 



n ^ ' 

Equation immediately suggests that, when building an azimuthal asymmetry of the type described in SecQ] if for 
some reason the leading term W°{p) gets cancelled, an asymmetry could arise from the projectile-nucleus interaction 
through the interference term G{p). 

Therefore, small contributions (p/ <^ p) can become important in an asymmetry. Consequently, we have to go back 
to the problem of the spin dependence of the projectile-nucleus interaction, which was initially neglected. 



III. SPIN-ORBIT EFFECTS 



The eikonal approximation is based on the assumption that a fast moving object experiences a roughly constant 
nuclear density but on the nuclear surface, where modifications of the projectile wave function can happen due to soft 
spin-orbit interactions. However, it is well known that the spin-orbit potential Vso is proportional to the gradient of 
the nuclear density p. For a projectile travelling along the usual z axis, the spin-orbit contribution ijjso to its wave 
function in the eikonal approximation reads 

i^soiz) - r dz' Vso{z') oc r dz' ^ = p{z) . (13) 

In other words, the modifications to the wave function induced by the spin-orbit potential do not depend on the 
details of the surface of the nucleus but, rather, on its bulk density. 

Since the projectile momentum is p — p z, it is natural to quantize the spin along the y direction such that the 
operator form of the spin-orbit potential is 

Vsoir)(T ■ r X p oc -^pxay , (14) 

where cr is the vector of 2x2 Pauli matrices. Consequently, the two spin components 'ip± of the incoming wave function 
(if the projectile has spin i) feel a spin-orbit potential of opposite sign; at large energies, Vso is also purely absorptive 
so that the net effect is that different polarization components of the incoming hadron beam feel different dampings. 
In the given conditions, the asymmetry should be reasonably independent from z, because it is created when the 
projectile enters the nuclear surface, and it is along the x direction: the nucleus is divided in two emispheres (e.g., 
a; > and a; < in the frame where the target is at rest and centered in the origin), each one being crossed by a flux 
of projectiles with opposite polarization. 

Since the spin-orbit damping adds to the nuclear damping described in the previous section, we can formally write 
the whole wave function as 



V'±(r) = ij{z){l±a{p)px) , 



(15) 



where il^{z) is defined in Eq. Q. The cocfhcient a{p) is largely energy dependent [s^. For sake of simplicity (and 
noting that the overall spin-orbit effect is not an increasing function of p at large energies |34j ^ . we rewrite the 
asymmetric factor as a{p)p = ?7(p). Since 7]R ^ 1 at large energies, we may approximate 1 ± 77a; « exp(±rya;) and 
write the components of the projectile wave function as 



V^±(r) 



Jpz-pi z±r]x 



(16) 



which correspond to two plane waves with complex momenta P± = P z ^ P^x — {p + ipi) z^irjx, respectively. 
We can now repeat the previous steps leading to Eq. ^ and write the scattering amplitude as 

a±{p + ipi,iri) = J dkA{{p + ipi)z T irix,k) . (17) 

Again, following Eqs. H9I12(I we can perform a Taylor expansion including only terms linear in pi and i], namely 

da±n 



a±niP ^ p + ipi,Px ^iri) ~ a±n{p) + ipi ■ 



dP 



. da±n 



pj=ri=0 



(18) 



which leads to the response 



where 



W±ip,pi,Ti) = W^{p) - 2pjlm [G^ip)] ± 2r,Im [G 



da 



±n 



'■±71 ' 



dP 



p/=r)=0 



, Gi{p)=j:[- 



'in 



da±, 



pj=rj=0 



(19) 



(20) 



IV. AZIMUTHAL ASYMMETRY IN UNPOLARIZED DRELL-YAN 



When building an azimuthal asymmetry from a hard event, the response W± must depend also on a variable, say 
^, which is characteristic of the asymmetric distribution of the final products. In the example of Fig. ^ ^ can be the 
angle (j) which identifies the azimuthal direction of the lepton pair production in the CoUins-Soper frame. By making 
this dependence explicit in W± , a typical feature of these asymmetric distributions for a free proton target (or when 
no distortion is taken into account, i.e. for pi = -q = 0, hence the superscript ° in the following formula) is that for 
each ^ there exists a correlated ^' such that 



M^^(^,p)-W^-(r,p) 



(21) 



In the previous example, the correspondence ^ — » ^' could be — > —(j) or — > + tt. In other words, there are always 
two alternative ways for defining and/or calculating an azimuthal spin asymmetry when the effect of nuclear medium 
is not included: fix ^ and reverse the spin of the projectile, or fix its spin and replace ^ by 

Then, it is natural to define the spin asymmetry for a polarized Drell-Yan event on a free proton target as 



Wl{^,p) - Wl{i,p) ^ Wl{^,p) - Wl{£,',p) 
Wl{i,p) + W°{(,,p) Wl{^,p) + Wl{^\p) 

In addition, we also define the "polarization averaged" asymmetry 



(22) 



A°^\ iwiitp) + W^f (e,p) ~ Wl{S,\p) - Wl{£,',p) 



(23) 



Because of Eq. (I21|l . the latter A° turns out to vanish. It means that those scattering amplitudes that contribute to 
Att simply average to zero. Effects of a different origin may create a nonvanishing A° . In fact, because of the nuclear 
damping summarized in Eq. H19|l , the analogue of Eq. (|21|l for a real nuclear target becomes 



(C, p, p/ , v) - iC ,P,Pi,rj) = 2pi Im [G^ (C , p) - G^ (C, p)] + 2,7 Im [G^ (C, p) + G^ (e' , p)] , 



(24) 



where the " undistorted" contributions cancel just because of Eq. (|21|l . If p/ and/or rj are not vanishing, then also 
Eq. 124|) is not vanishing. Consequently, the generalization of A° to the case of a real nuclear target reads 



A 



W+{^,p,Pi,r]) + W-{^,p,pi,r]) - W+{^',p,pi,Tj) - W-{^',p,pi,?]) 



- I 2p,Im [G^(C',p) + G^(e',p) - G^(C,p) - G^(e,p)] 



+2,7 Im [Glitp) - Gl{i,p) - Gl{i',p) + G^(C',p)] 



(25) 



When nuclear distortion is neglected or a free proton target is considered, i.e. for p/ = = 0, then Eq. H21(l holds 
and A = A° = Q. It means that for a nuclear target the projectile can equivalently follow different paths (^ and ^') 
inside the nucleus. When including the nuclear distortion, from Eq. H25II we get in general that A^Q: because of the 
nuclear damping, the two paths ^ and ^' are not equivalent and the nucleus acts as a polarizer for the projectile. 

We can now put this picture in the physical context described in SecQ] If the Lam- Tung sum rule were correct, the 
cross section would not depend on the azimuthal (j) distribution of the final lepton pair in an unpolarized Drell-Yan 
event in nuclear targets. This is somewhat the same content of Eq. (|21|l and of the condition A° = for ^ = (j). 



However, as already stressed in Sec. the Lam- Tung sum rule is valid only for the direction of annihilation along 
the z axis (see Fig.^l, i.e. for collinear partons in the hard event. Violation of the sum rule emphasizes the role of 
transverse dynamics of partons inside hadrons. A similar azimuthal asymmetry can be obtained by considering that 
the transverse momentum of the Drell-Yan pair originates from the transverse momentum of the colliding hadron 
when it is deflected by nuclear damping (particularly by the "transversely asymmetric" spin-orbit effect). This sort of 
duality that is intriguingly established between the partonic and hadronic interpretations of azimuthal asymmetries in 
terms of transverse (spin) degrees of freedom, might suggest that spin-orbit effects play a role also at the elementary 
level in characterizing the transverse dynamics of partons inside hadrons. 

In conclusion, when a Drell-Yan hard event is considered in a collision of a hadronic unpolarized projectile with an 
unpolarized nuclear target, the nuclear damping and the soft spin-orbit interaction happening on the nuclear surface 
can distort the projectile wave function and produce an azimuthal asymmetric distribution of the final lepton pair 
(v4 7^ 0) even if there is no asymmetry in the corresponding collision on a free proton target (A° — 0). The necessary 
condition for this to happen is that for the corresponding Drell-Yan with a polarized projectile and/or proton target 
the transverse spin asymmetry is not vanishing (Att ^ 0), because both observables are expressed in terms of the 
same scattering amplitudes a„. 



V. DISCUSSION AND PERSPECTIVES 



Given the results summarized in the last paragraph of previous section, the question arises about the interpretation 
of the experimental data of Rcfs. and of possible future Drell-Yan experiments with (anti)proton beams 

at the foreseen HESR at GSI: when a nuclear target is employed, how big is the distortion of the beam wave function 
at a given kinematics? And how can this effect be disentangled from the genuine azimuthal asymmetry due to the 
hard event? 

Concerning the first question, the main problem is the lack of knowledge about the microscopic mechanisms that 
build up the elementary scattering amplitudes a„ in Eq. . We can attempt some educated guess based on simpli- 
fying hypothesis and on the general properties of known parton distributions and of hadron-hadron phenomenology. 

For sake of simplicity, let W° in Eq. jnj be dominated by a single amplitude a(P, Px) = -^ia^ with daJdPi^^) ~ 
and « at p/ = 77 = 0. Then, 



W^°(p) 



(26) 



where i — L,T correspond to longitudinal and transverse components with Pl — P = p 
We define the ratios 



ipi and Pt — Px = irj. 



ri 



W°{p) 



W°{p) 



(27) 



that should tell us about the relative size between the asymmetry of final lepton pairs generated by the collision of an 
unpolarized hadron beam on a nucleus, and by the collision of the same beam, but polarized, on a free proton target. 
In other words, r/ and r^o should tell us how big is the asymmetry produced by nuclear damping and spin-orbit, 
respectively, with respect to the asymmetry produced by the genuine polarized Drell-Yan hard event. 

In order to estimate r/ and r^o, we observe that the parton densities significantly change over the scale of their 
parent hadron size Rh ^ 0.2 0.5 fm (which corresponds to a momentum scale of p;i ~ 1 0.4 GeV/c, respectively). 
We assume that also the amplitude a changes significantly only over the same scale; neglecting fine details and looking 
only to the order of magnitudes, we can roughly conclude that 



G\P) 
W°{p) 



1 da^^ 



pi=ri=0 



— Rfi 

Ph 



(28) 



In other words, we safely assume that the relevant mechanisms involve correlations at most over the scale Rh excluding 
long-range interactions that would imply large derivatives of a a,t pi = j] — and, consequently, large values of G". 
In this hypothesis, the ratios of Eq. l|77|l become r/ ^ 2pi/ph and Tso ~ 2r\lph., respectively. As already mentioned 
in Sec.m experimental data about nuclear transparency on ^^C seem to suggest pi ^ 50 MeV/c. This value can be 
interpreted as the mean free path of a hadron travelling in the nuclear medium, and it should not depend on the mass 
number but for very light nuclei. The corresponding ratio r/, then, could be in the range 5 10 %. 

In order to estimate 77, it is useful to consider the nuclear analyzing power Ajsi, i.e. the asymmetry between cross 
sections at a given transferred momentum for spin up and spin down projectiles colliding on an unpolarized nucleus. 



If the typical range associated with spin-orbit effects is given by the surface thickness, we can approximate this scale 
in a nucleus again with Rh- Hence, we can assume that (An) « t] ^ T]/ph, where (A^) is the analyzing power 
averaged over the transferred momenta relevant to the ISI occurring during the propagation of the projectile inside 
the nuclear target, typically in the range 0.5 GeV/c. The problem of estimating r^o is, therefore, translated into 
the determination of the typical size of (A]\j). In the literature, recent measurements and attempts of theoretical 
interpretations have been published for elastic scatterings at energies beyond 10 GcV with proton-proton |39t l4rl . 
proton-^^C .41, .42,|0|, and antiproton-proton systems (for a collection of older experimental data, see Ref. '42j 
and references therein). Results for An with proton and ^^C targets are similar in magnitude and shape and they show 
a steep decrease from 10% at 10 GeV down to 2-^4% at 20 GeV, with somewhat higher values at larger energies around 
30 GeV but with l arge error bars. At large energies, the antiproton-proton results have also a similar magnitude to 
the previous ones |39l |. However, extrapolation of these findings to other nuclear targets (maybe, isoscalar nuclei 
like ^^C) should be taken with great care, since the reproduction of "nuclear" An in terms of proton-proton and 
antiproton-proton amplitudes is not theoretically well established For transferred momenta in the range of 

interest 0.5 GcV/c, A^ already gets its maximum value for almost all beam energies and the other values are not 
much smaller than that. At larger angles. An can reach large values at all energies, but this situation is not relevant 
here and it won't be considered. In conclusion, with the above caveat we can assume (An) ~ An and finally deduce 
that the approximate size of the spin-orbit nuclear damping is in the range 4 10%, depending on the projectile 
energy. 

The previous estimates about r/ and r^o are rather conservative, since they are based on hard interactions occurring 
at most within the short range Rh- However, we cannot exclude the case of a scattering amplitude a having a much 
steeper dependence on small values of P{Px) because of correlations at a longer range. For example, if a has a 
resonance-like behaviour with a pole as 



a(fc) = —— , (29) 

K — 27 



where k = P or Px, then from Eq. H26(l we get 



_ 2pjGHp) _ 2pz _ 2r,G^{p) _ 2^ 

W"ip) 7 ' wip) ~ 7 ■ ^ ' 

In the "short-range" hypothesis, 7 ^ 1/Rh- But if a represents an intermediate state with a narrow width, then the 
nuclear damping is largely enhanced. This is not surprising, since a narrow-width resonance would propagate inside 
the nucleus with a mean free path '^1/7 increasing the importance of nuclear effects. A possible example of this 
phenomenon is represented by the Drell-Yan process at the J/^ mass [26j, but involving a nuclear target. Since this 
resonance is very narrow and it is present in all possible amplitudes, if the process is studied via a hadron-nucleus 
collision we would have the standard conditions of nuclear shadowing, i.e. the propagation of a quasi-real particle 
through large distances inside the nucleus. 

In summary, nuclear damping effects in hadron-nucleus collisions will affect azimuthal asymmetries from Drell-Yan 
hard events most likely by few percents, unless quasi-real intermediate states are involved in the propagation of the 
hadron projectile inside the target nucleus before the hard event. This condition is not so unlikely, if higher-twist 
contributions are relevant in Drell-Yan processes. Indeed, if Mp/M is related to the ratio between subleading and 
leading twists (with Mp the proton mass), a simulation in the typical kinematics foreseen at HESR at GSI [l^l 
reveals that the cross section has a steep M dependence with most events concentrated at the lowest allowed M = 4 
GeV, where Mp/M « 25%. However, new measurements of Drell-Yan events are of course needed to confirm and 
substantiate previous arguments, since the available data pertain a completely different regime j^lfll. [in|.lllj. where, 
incidentally, large error bars in the azimuthal asymmetries do not exclude nuclear effects when changing from the 
deuteron to the tungsten targets (see Fig. 8 of Ref. 0). 

Going back to the beginning of the section, we consider the second question about the possibility of disentangling 
the effect of nuclear damping from the genuine azimuthal asymmetry produced by a Drell-Yan hard event. Since 
spin-orbit effects are mostly associated with a single quasi-elastic projectile-nucleon scattering, detection of the Drell- 
Yan lepton pair in coincidence with a recoil proton with small longitudinal momentum (< 300 MeV/c) and large 
transverse momentum (400 800 MeV/c) would allow for a quantitative check of the importance of ISI (responsible 
for this large-angle emission) and for an indirect determination of the polarization of the projectile before the hard 
event, because this would determine the direction of the recoil proton momentum. With this tagging technique, it 
should be possible, in principle, to study a polarized Drell-Yan event from unpolarized beams and targets. 



VI. CONCLUSIONS 



Since several years, a series of measurements of Drell-Yan events in unpolarized high-energy hadronic coUisions on 
nuclei have been published in the literature 0, 0, that still await for a coherent quantitative explanation. In 
fact, the differential cross section shows an unexpected largely asymmetric azimuthal distribution of the final Icpton 
pair with respect to the production plane (see Fig. At present, the most promising interpretation identifies the 
source of the observed cos 2<j) asymmetry in a leading-twist contribution to the cross section containing a peculiar 
distribution function, that describes the intrinsic transverse momentum distribution of transversely polarized partons 
inside unpolarized hadrons ^13] . 

In this paper, we have explored the consequences of the fact that the above experiments involve target nuclei, namely 
that the hadronic beam (pions or antiprotons) has to cross the nuclear surface and travel inside the nuclear matter 
for some length (related to its energy) before colliding on a bound nucleon to produce the Drell-Yan hard event. The 
wave function of the projectile is distorted by the nuclear damping and by spin-orbit interactions happening on the 
nuclear surface, producing also an azimuthal asymmetry in the distribution of the Drell-Yan lepton pair; this effect 
should be added to the one originating from the elementary hard event. The necessary condition for this to happen 
is that for the corresponding Drell-Yan with polarized projectile and/or free proton target the asymmetry produced 
by flipping one of the spins is not vanishing, since both asymmetries are expressed in terms of the same elementary 
scattering amplitudes. 

Using some educated guess based on simplifying hypothesis and on the general properties of known parton distri- 
butions and of hadron-hadron phenomenology, we estimated that the nuclear damping should affect the azimuthal 
asymmetry from Drell-Yan hard events most likely by few percents, unless quasi-real intermediate states are involved 
in the propagation of the hadron projectile inside the target nucleus before the hard event. This condition is not so 
unlikely if nuclear targets will be possibly adopted at the HESR ring at GSI, where in the considered kinematics either 
higher-twist contributions are relevant or the invariant mass of hadronic resonances with narrow width is reached in 
the annihilation (see Refs. [H IH [H Hi IS El 113 for further details). 

It should be possible to isolate " nuclear" asymmetries from " hard" asymmetries by considering a Drell-Yan process 
where the final lepton pair is detected in coincidence with a recoil proton with small longitudinal momentum and large 
transverse momentum. In fact, only the rescatterings of the hadron beam propagating inside the nuclear medium 
would be responsible of such large-angle emission. It is interesting to note that, as the cos 20 asymmetry could be 
produced by the transverse dynamics of partons inside hadrons leading to a violation of the Lam- Tung sum rule 
similarly the same asymmetry can be obtained by considering the transverse momentum of the colliding hadron when 
it is deflected by nuclear damping (particularly by the "transversely asymmetric" spin-orbit effect). This parallelism 
might suggest that spin-orbit effects play a role also at the elementary level in characterizing the transverse dynamics 
of partons inside hadrons. 



[1] G. Bunce et al., Phys. Rev. Lett. 36, 1113 (1976). 

[2] D. L. Adams et al. (FNAL-E704), Phys. Lett. B264, 462 (1991). 

[3] A. Bravar (Spin Muon), Nucl. Phys. Proc. Suppl. 79, 520 (1999). 

[4] J. Adams et al. (STAR), Phys. Rev. Lett. 92, 171801 (2004), hep-ex/0310058. 

[5] M. Radici and G. van der Steenhoven, CERN Courier 44, 51 (2004). 

[6] R. L. Jaffe (1996), proceedings of the Ettore Majorana International School on the Spin Structure of the Nucleon, Erice, 
Italy, 3-10 Aug 1995., hep-ph/9602236. 

[7] V. Barone and P. G. Ratcliffe, Transverse Spin Physics (World Scientific, River Edge, USA, 2003). 

[8] S. Falciano et al. (NAIO), Z. Phys. C31, 513 (1986). 

[9] M. GuanziroU et al. (NAIO), Z. Phys. C37, 545 (1988). 
[10] J. S. Conway et al., Phys. Rev. D39, 92 (1989). 
[11] E. Anassontzis et al., Phys. Rev. D38, 1377 (1988). 
[12] J. C. Collins and D. E. Soper, Phys. Rev. D16, 2219 (1977). 
[13] D. Boer, Phys. Rev. D60, 014012 (1999), hep-ph/9902255. 
[14] A. Brandenburg, O. Nachtmann, and E. Mirkes, Z. Phys. C60, 697 (1993). 
[15] C. S. Lam and W.-K. Tung, Phys. Rev. D21, 2712 (1980). 

[16] A. Brandenburg, S. J. Brodsky, V. V. Khoze, and D. MuUer, Phys. Rev. Lett. 73, 939 (1994), hep-ph/9403361. 
[17] K. J. Eskola, P. Hoyer, M. Vanttinen, and R. Vogt, Phys. Lett. B333, 526 (1994), hep-ph/9404322. 
[18] E. L. Berger and S. J. Brodsky, Phys. Rev. Lett. 42, 940 (1979). 

[19] R. J. Fries, A. Schafer, E. Stein, and B. Muller, Nucl. Phys. B582, 537 (2000), hep-ph/0002074. 

[20] B. Z. Kopeliovich, J. Raufeisen, A. V. Tarasov, and M. B. Johnson, Phys. Rev. C67, 014903 (2003), hep-ph/0110221. 



[21] PANDA CoUab., Letter of Intent for the Proton-Antiproton Darmstadt Experiment, Spokesperson: U. Wiedner (2004), 

http://www.gsi.de/documents/DOC-2004-Jan-115-l.pdf. 
[22] PAX CoUab., Letter of Intent for Antiproton-Proton Scattering Experiments with Polarization, Spokespersons: F. Rath- 

mannn and P. Lenisa (2004), http://www.fz-juelich.de/ikp/pax/. 
[23] P. Lenisa ct al. (PAX) (2004), hop-cx/0412063. 

[24] ASSIA CoUab., Letter of Intent for A Study of Spin-dcpcndcnt Interactions with Antiprotons: the Structure of the Proton, 

Spokesperson: R. Bertini (2004), http://www.gsi.de/docuinonts/DOC-2004-Jan-152-l.ps. 
[25] A. V. Efremov, K. Goeke, and P. Schweitzer, Eur. Phys. J. C35, 207 (2004), hep-ph/0403124. 
[26] M. Anselmino, V. Baronc, A. Drago, and N. N. Nikolaev, Phys. Lett. B594, 97 (2004), hep-ph/0403114. 
[27] M. Radici (2004), hcp-ph/0407085. 

[28] N. C. R. Makins et al. (the NE18 Collaboration), Phys. Rev. Lett. 72, 1986 (1994). 

[29] R. J. Glauber, Lectures in Theoretical Physics, vol. 1 (Interscience Publ., New York, 1959). 

[30] G. D. Alkhazov, S. I. Belostotsky, and A. A. Vorobycv, Phys. Rep. 42, 89 (1978). 

[31] A. Bianconi and M. Radici, Phys. Lett. B363, 24 (1995). 

[32] A. Bianconi and M. Radici, Phys. Rev. CSS, R563 (1996). 

[33] A. Bianconi and M. Radici, Phys. Rev. C54, 3117 (1996). 

[34] S. Boffi, C. Giusti, F. D. Pacati, and M. Radici, Electromagnetic Response of Atomic Nuclei, vol. 20 of Oxford Studies in 

Nuclear Physics (Oxford University Press, Oxford, 1996). 
[35] A. Bianconi and M. Radici, Phys. Rev. C 56, 1002 (1997). 

[36] Y. S. Golubeva, L. A. Koudratyuk, A. Bianconi, S. Bofli, and M. Radici, Phys. Rev. C 57, 2618 (1998). 
[37] A. S. Davydov, Kvantovaja Mechanika (Nauka, Moscow, 1981). 
[38] J. P. Ralston and D. E. Soper, Nucl. Phys. B152, 109 (1979). 
[39] N. Akchurin et al. (E581/704), Phys. Lett. B229, 299 (1989). 

[40] N. H. Buttirnorc, B. Z. Kopcliovich, E. Leader, J. Soffer, and T. L. Trueman, Phys. Rev. D59, 114010 (1999), hep- 
ph/9901339. 

[41] J. Tojo et al., Phys. Rev. Lett. 89, 052302 (2002), hep-ex/0206057. 
[42] B. Z. Kopeliovich and T. L. Trueman, Phys. Rev. D64, 034004 (2001), hep-ph/0012091. 
[43] O. V. Selyugin and J. R. Cudell, Czech. J. Phys. 54, B161 (2004), hep-ph/0301048. 
[44] A. Bianconi and M. Radici (2004), hep-ph/0412368. 



